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The paper proposes a new model of spin dynamics which can be treated as a model of 
sociological coupling between individuals. Our approach takes into account two different 
human features: gregariousness and individuality. We will show how they affect a psycho- 
logical distance between individuals and how the distance changes the opinion formation 
in a social group. Apart from its sociological aplications the model displays the variety 
of other interesting phenomena like self-organizing ferromagnetic state or a second order 
phase transition and can be studied from different points of view, e.g. as a model of 
ferromagnetic fluid, complex evolving network or multiplicative random process. 
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1. Introduction 

Interdisciplinary research has been drawing much attention in the last decades. Mod- 
els and methods developed in theoretical physics proved to be fruitful in studying 
complex systems composed of relatively simple mutually interacting elements 
and coming from domains as diverge as neural networks ^ , disease spreading pop- 
ulation dynamics ^, etc. But the range of the investigations goes also beyond the 
natural sciences and includes problems from sociology or economy, like pedestrian 
motion and traffic ^, migrations '''^ or financial crashes ®. Another important sub- 
ject of this kind is the process of opinion formation in social groups. One way of 
its quantitative description consists in a macroscopic approach based on the master 
equation or the Boltzmann-like equations for global variables Alternatively, 
by making some sociologically motivated assumptions on the mechanisms of inter- 
actions between individuals " microscopic" models are constructed and investigated 
numerically or analytically by means of methods known from statistical physics 
I2,i3_ Qj^g concludes that the variety of the emerging physical collective phenomena 
has much in common with the complex social processes. 

In particular, Nowak, Szamrej and Latane created a simple model based upon 
the successful theory of social impact in human societes first introduced by La- 
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tane in 1981 In the simplest form their model characterizes the strength of the 
psychological coupling between the individuals by two qualities: persuasiveness and 
support. The former describes the ability of one individual to persuade the other 

one to change his/her opinion. The later describes the ability of one individual to 
support the other one in his/her opinion. Different variants of the model were ex- 
plored numerically 15,16,17^ ^^^^ many of the observations were than explained in 
the framework of a mean field approach 18.19,20 ^^^^ ^j^g Landau theory 

Here wc would like to present a rather different approach to describe psycho- 
logical coupling. Instead of persuasiveness and support we will study the effect of 
gregariousness and individuality. We will show how the two features may affect 
psychological distance between individuals and how the distance changes opinion 
formation in the society. Finally, we will show that our model could be mapped to 
ferromagnetic fluid not in Euclidean but in a social space. 

2. The model 

Our system consists of N individuals (members of a social group); we assume that 
each of them can share one of two opposite opinions on a certain subject, denoted 
as ai = ±1, i = 1, 2, ...N. The Hamiltonian of the model reads: 



Individuals can influence each other with the strength Jij{t) which can be un- 
derstood as an inverse of their distance in a social space. The above means that a 
stronger impact corresponds to a shorter distance. We assume that social distances 
are changing in time and we put on the following dynamics of the strength Jij{t): 



The parameter rj > is responsible for continous growth of the social strength 
and can be identified as gregariousness of i-th individual which leads to tightening 
of ties with other people. In other words, people from their nature seek the company 
of others. The parameter a > describes another natural tendency of people which 
is a need to be different than a surrounding crowd, i.e. it reflects the inclination of 
an individual to demonstrate his/her individuality. 

For completeness of the model we assume as an initial condition any positive val- 
ues of Jij{t = 0). The condition assures that during the system evolution couplings 
are always positive in the most interesting range of parameters rj and a. 

Now, let us concentrate on the phase diagram for the presented model (fig. 1). It 
is divided into four different regions by three curves. The curve 1 is the most obvious 
one. For every set of parameters above this curve, i.e. for rj > a, coupling strengthes 
will increase to infinity in exponential way. Parameters 77 < a — 1, limited by the 
curve 3, also make the system unstable, but now coupling strengthes can become 
negative. It means that in every step J will change its sign and | J| will diverge to 
infinity. The stable region lies between two curves 1 and 3. 




(1) 



Jij{t + 1) = Ji,j{t){l + 7]- OLCFiaj). 



(2) 
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Fig. 1. Phase diagram of the model (l)-(2). F - ferromagnetic phase, P - paramagnetic phase. 
Detailed explanation in the text. 

To explain the curve 2 let us concentrate on a single coupling and investigate 
the following process: J{t) — » J{t + 1) — » J{t + 2) = J{t). As one can see from eq. 
2, in every step J grows or decreases by some well defined value. The above process 
is the simplest one in which J stays at some fixed level, i.e. it grows and then it 
decreases to the same value. Let us assume that J{t + 1) > J{t) (the opposite case 
is analogous). It means that J{t+1) = J{t){l+r] + a). Then, in next step J should 
decrease ,i.e. J{t + 2) = J{t+l){l + r] — a). Prom above one can obtain the following 
condition: 

r?e = VT+^ - 1. (3) 

It is easy to see that a corresponding critical condition for a generalized process 
J{t) — *■ J{t + 1) J(t + 2) ^ ... —>■ J(t + 2n) = J{t) is equivalent to the condition 
(3). 

If we set 7] below the critical value rjc then the above process will result in 

decreasing J to zero. Then it is obvious that two spins for which this coupling J 
is investigated become disconnected and independent, what leads to paramagentic 
state. 

To complete discussion of the phase diagram one has to note that it is useful to 
assure a <C 1. If one rewrites the eq. (2) in the following form: 

Jij{t+1)-J,j{t) ^ A A, it) ^ri_ 

aJi,j{t) J^,J{t)aAt a ''''''' ^ ' 

then one can see that a plays the role of time scale. It means that for large a the 
succcsive values of J are very distant (J changes very fast) and spin dynamics can 
not follow to compensate changes of J. It manifests itself in long time observed 
paramagnetic states interrupted by long time observed ferromagnetic states. 

In summary, the interesting from the sociological point of view range of param- 
eters is ricritical < r? < a < 1. 
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The dynamics of changes of individual's opinion is given by a simple Monte 
Carlo procedure based on the Metropolis algorithm. A temperature T given in 
the algorithm may be interpreted as a "social temperature" describing degree of 
randomness in the behavior of individuals, but also their average volatility. The 
procedure consists of two steps. In the first step we update states of N randomly 
chosen individuals. In the second step we update coupling strengths for all nodes 
according to cq. (2). 

As we will show for a wide range of parameters rj and a, regardless of choosing 
a temperature the system tends to be in a ferromagnetic regime. It means that 
despite a tendency to manifest individuality most of individuals interact with the 
other people who share the same opinion. 

3. Results 

A typical dependence of magnetization per spin \m\ on system parameters rj/a is 
shown at fig. 2. Considering r] as an order parameter, continous (second order) phase 
transition occurs for r]c given by eq. (3). Open points presented at fig. 1 obtained 
from simulations confirm that the above derivation is correct. 

One can see from the fig. 4 that the absolute value of the mean magnetization 
is an (increasing) function of ri/a but it is completely independent on the system 
temperature T. The fact can be understood as follows. According to cq. (4) we can 
write the following equation for the mean value for the logarithm of Jij 

1 /dlnJij\ V , \ 

. 8 
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Fig. 2. Dependence of average magnetization per spin on the system parameter rj/a. 

However, the mean value of the product {(JiCTj) is related to the mean system 
magnetization (m) which on the other hand is a certain function of the Boltzmann 
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factors exp(— Jj_j/r). Thus taking andvantage of the mean field approximation we 
can write that 

When the system is close to equilibrium the left hand side of the last equation in 
average equals to zero and the equation simplifies to the following relation 

(J.,,)/r = 5-i(ry/a). (7) 

In this sense the average value of the coupling constant is always proportional to 
the temperature and a function of a ratio rj/a. 

Since m is a function of the ratio {Jij)/T thus it only depends on the ratio r]/a 
and does not depend on the system temperature. The numerical confirmation of the 
statement is presented at fig. 3. 
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Fig. 3. Average coupling strength in the system as a function of temperature. Log-log fit is 
presented by straight lines with slopes 1.03 for open points and 1.01 for filled jjoints. Inset: the 
same dependence in linear scale. Slopes of two fitted linear functions represent a pioportionaUty 
fa<;tor given by function g~^(r)/a) from eq. (7). 

To illustrate the effect one can perform the following procedure. First, set values 
of the parameters rj and a, and also temperature Ti of the system. (These values re- 
sult in a given spin magnetization depicted in fig. 4 by m* , and they also correspond 
to a certain distribution of couplings between spins P{ J).) Now we would like to 
reconstruct the whole magnetization curve passing through the point (Ti,to*). Un- 
fortunately, as we have already shown any change of T modifies coupling strengthes, 
which in consequence modify the shape of the magnetization curve. Reconstruction 
of m{T) will be possible if we freeze P{J), i.e. we make each Jij constant. 

Once we have determined m{T), we restore dependence of J on T and change 
temperature to new value T2 > Tci- A new pair (T2,m*) determines a new critical 
temperature Tc2- The curve at the fig. 4a adjusts to new conditions and transforms 
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a) b) 




Fig. 4. Schematic picture of processes occuring in the system during the changing of the tem- 
perature T and rj/a = const; a) equilibrium (initial) state of the system; b) The change of the 
temperature T from Ti to T2 forces magnetization curve to reshape to new equilibrium conditions. 

to the shape shown at the fig. 4b (obtained by the same method as before). It 
means that regardless of a choice of the system temperature we are always below 
the critical temperature, i.e. in the ferromagnetic state. 

The sociological conclusion could be as follows: regardless of a "social temper- 
ature" people always try to correlate their opinions with others (create groups of 
interest). This tendency to share the same opinion with other people, regardless of 
some external forces, make us, people, so resistant to trials of despots to make the 
people unorganized and disoriented. Of course the parameter rj/a characterizes our 
own (not social) point of view which gives us some independency respecting other 
people opinion. 

From the point of view of complex networks domain 22,23 j-j- jg interesting to 
consider the model as a weighted network, where nodes correspond to individuals 
and links have assigned weights equal to coupling strength. One of nontrivial obser- 
vations is a distribution of coupling strengthes P{J) which is presented at fig. 5. As 
one can see for large temperature T the distribution has a form of power law with 
the exponent 7 w 0.85. 

It seems that there should be a strong relation between the observed power-law 
distributions and distributions obtained due to a more general class of multiplicative 
random processes 24,25_ j£ ia,ct, one can easily find some similarity of eq. (2) to eq. 
(1) in The differences occur when one takes into account the temperature and its 
infiuence on distributions at fig. (5). We suspect that the model studied by us settles 
somewhere between two multiplicative random processes studied in and . This 
hypothesis is still under investigation and the results will be published elsewhere. 

Now let us draw attention to similarity of the presented model to magnetic 
fluids which are widely studied for the last thirty years Magnetic fluids are 
described by interacting molecules with both translational and spin degrees of free- 
dom. They interact due to weak long-ranged exchange interactions in addition to 
spin-independent isotropic attractive forces. The most simple physical parameter 
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Fig. 5. Distribution of coupling strengths for different values of parameters rj/a and T. 



used in phase diagrams of magnetic fluids has a form 

/ (t>ex{r)dr 



R = 



(8) 



where 4)ex (^) describes exchange integral, 4>attr {t^ some attractive-type integral 
and r is Euclidean distance between molecules. One can easily see a correspondence 
between the above parameter and the main parameter of our model rj/a. In that 
sense 77 and a have respectively attractive and spin-dependent properties, and our 
distance in social space reflects Euclidean distance r. 



4. Conclusions 

In this paper we propose a new model of spin dynamics which could be treated 
as a model of sociological coupling between individuals. Apart from its sociological 
aplications the model displays the variety of other interesting phenomena like self- 
organizing ferromagnetic state or a second order phase transition and can be studied 
from different points of view, for example as a model of ferromagnetic fluid, complex 
evolving network or multiplicative random process. 
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